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V-ry  ^  ABSTRACT 

/Ai5  uiicikjhtr  “ 

~Wer  oonsider^he  flow  of  two  immiscible  fluids  lyinq  between  concentric 
cylinders  when  the  outer  cylinder  is  fixed  and  the  inner  one  rotates.  The 
interface  is  assumed  to  be  concentric  with  the  cylinders  and  gravitational 

TU 

effects  are  neglected.  -We- present  a  numerical  study  of  the  effect  of 
different  viscosities,  different  densities  and  surface  tension  on  the  linear 

flu  r 

stability  of  the  Couette  flow,  pair’  results  indicate  that  with  surface 
tension,  a  thin  layer  of  the  less  viscous  fluid  next  to  either  cylinder  is 
linearly  stable  and  that  it  is  possible  to  have  stability  with  the  less  dense 
fluid  lying  outside.  The  stable  configuration  with  the  less  viscous  fluid 
next  to  the  inner  cylinder  is  more  stable  than  the  one  with  the  less  viscous 
fluid  next  to  the  outer  cylinder.  The  onset  of  Taylor  instability  for  one- 
fluid  flow  may  be  delayed  by  the  addition  of  a  thin  layer  of  less  viscous 
fluid  on  the  inner  wall  and  promoted  by  a  layer  of  more  viscous  fluid  on  the 
inner  wall. 
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SIGNIFICANCE  AND  EXPLANATION 


Steady  two-fluid  flow*  of  imiscible  density-matched  fluids  with 
different  viscosities  arise  in  applications  such  as  the  pipeline  transport  of 
oil  with  the  addition  of  water,  the  formation  of  bicomponent  fibers  such  as 
nylons  and  modelling  of  the  Earth's  mantle.  Such  flows  are  typically 
nonunique,  even  when  the  speeds  involved  are  slow.  However,  experiments 
usually  result  in  very  stable  unique  arrangements.  In  order  to  get  an  idea  of 
which  interface  positions  are  allowed  for  the  'Taylor  Problem'  we  study  the 
linear  stability  of  the  concentric  arrangement.  We  find  that,  in  the  absence 
of  gravity,  a  thin  lubrication  layer  of  the  less  viscous  fluid,  lying  next  to 
either  cylinder,  is  linearly  stable.  This  aqrees  with  the  experimental 
observations.  The  fact  that  the  less  viscous  liquid  tends  to  shield  the  more 

viscous  fluid  from  shearing  suggests  that  the  use  of  two  such  fluids  in 

lubrication  may  be  economical.  We  also  study  numerically  the  effect  of 

different  densities  and  surface  tension  and  find  that  the  stabilizing  effect 
of  this  viscosity  stratification  can  even  overcome  a  destabilizing  density 
difference:  the  arrangement  with  the  less  dense  fluid  outside  can  be  stable 
if  it  is  also  the  less  viscous  fluid  and  if  this  outer  layer  is  thin. 
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TWO-FLUID  COUETTE  FLOW  BETWEEN  CONCENTRIC  CYLINDERS 
Yuriko  Renardy  and  Daniel  D.  Joseph* 

Introduction 

We  consider  linear  stability  of  the  flow  of  two  immiscible  fluids  separated  by  an 
interface,  lying  between  concentric  rotating  cylinders.  In  each  fluid,  the  Navier-Stokes 
equations  for  steady  flow  are  assumed  to  hold.  If  we  prescribe  the  ratio  of  the  total 
voIum  occupied  by  each  fluid,  then  the  interface  is  an  unknown,  across  which  the  velocity 
and  normal  and  shear  stresses  are  to  be  continuous.  If  the  fluids  have  equal  or  nearly 
equal  densities,  then  a  continuum  of  interface  positions  are  allowed  (Joseph,  Renardy  and 
Renardy,  1983)*  However,  this  non-uniqueness  is  not  borne  out  by  the  experiments  of 
Joseph,  Nguyen  and  Beavers  (1983)  who  use  water  and  various  oils  as  the  two  fluids  in  an 
apparatus  with  the  outer  cylinder  fixed.  When  the  inner  cylinder  is  rotated  at  even 
moderate  speeds,  gravity  effects  appear  negligible  and  a  pattern  consisting  of  two  types 
of  cells  is  usually  observed.  One  type  consists  mostly  of  oil  rollers  stuck  to  the  inner 
cylinder  and  rotating  almost  like  a  solid  body,  lubricated  by  a  thin  layer  of  water  at  the 
outer  cylinder.  The  second  type  consists  mainly  of  water  cells  undergoing  Taylor  vortex 
motions.  These  cells  extend  from  the  inner  to  the  outer  cylinder  but,  in  some 
experiments,  are  covered  by  a  thin  layer  of  oil  at  the  outer  cylinder.  The  two  types  of 
cells  alternate  along  the  length  of  the  cylinder.  This  flow  is  one  of  many  steady 
bicomponent  flows  where  a  study  of  the  selection  mechanism  for  the  arrangement  of  the 
fluids  must  be  made.  One  way  to  study  selection  is  to  study  stability  and  in  this  paper 
we  study  stability  by  computing  eigenvalues  for  the  spectral  problem  associated  with  the 
linear  theory. 

The  equations  for  our  numerical  computations  are  given  in  Part  I.  Some  asymptotic 
results  for  short  waves  are  presented  in  Part  II  following  the  ideas  of  Hooper  and  Boyd 
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( 1993).  They  consider  unbounded  Couette  flow  but  their  Method  of  enelyeie  applies  locally 
et  any  interface  with  a  viscosity  jump.  Hoopsr  and  Boyd  showed  that  in  the  absence  of 
surface  tension,  the  flow  is  unstable  to  sufficiently  short  waves  which  have  wave-vectors 
parallel  to  the  basic  flow.  The  growth  rates  of  these  disturbances,  however,  tend  to  sere 
as  the  waves  get  shorter.  They  found  that  these  short-wave  instabilities  are  not 
suppressed  by  viscosity  as  they  are  in  one-fluid  flows,  but  by  surface  tension.  A  density 
difference  can  stabilise  or  destabilise  then  but  not  as  effectively  as  surface  tension. 

In  Part  II  t  we  have  given  a  similar  analysis  for  disturbances  whose  wave-vectors  are 
perpendicular  to  the  basic  flow.  Our  results  Include  surface  tension,  density  differences 
and  centrifugal  effects.  We  have  correlated  our  numerical  results  with  the  asynptotlc 
females. 

When  surface  tension  is  effective,  the  longer  waves  can  cause  instability  and  if 
periodic  boundary  conditions  are  imposed,  then  this  yields  a  familiar  type  of  instability 
in  which  the  interaction  of  a  finite  number  of  modes  determines  what  type  of  solutions 
bifurcate  from  the  unstable  one  (M.  Renardy  and  Joseph,  1983).  However,  when  surface 
tension  is  not  effective,  then  we  have  an  unusual  instability  in  which  the  flow  is 
unstable  to  all  short  waves  below  a  certain  critical  else.  This  type  of  Instability  may 
play  a  role  in  the  formation  of  emulsions. 

In  Part  III,  we  give  numerical  results  for  two  situations i  low  Reynolds  numbers  and 
Taylor  numbers  near  a  critical  one.  We  find  that  in  the  preferred  configurations,  a  thin 
layer  of  the  less  viscous  fluid  may  lie  next  to  either  cylinder.  Our  results  contradict 
the  selection  principle  based  on  minimising  the  viscous  dissipation  in  the  restricted 
class  of  annular  layers  of  two  fluids  which  do  not  vary  along  the  axis  of  the  cylinder. 

The  solution  of  this  minimisation  problem  (Joseph,  Nguyen  and  Beavers,  1983)  consists  of 
the  less  viscous  fluid  lying  on  the  inner  cylinder,  no  matter  which  cylinder  rotates.  In 
fact,  our  numerical  results  indicate  that  a  narrow  stable  layer  of  the  less  viscous  fluid 
on  the  inner  cylinder  is  "more”  stable  than  that  on  the  outer  cylinder.  It  is  also  of 
Interest  that  the  instability  leading  to  Taylor  vortices  in  one-fluid  flows  may  be 
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nullified  by  addin?  a  lubrication  layar  of  the  less  viscous  fluid  at  the  inner  cylinder 
and  that  instability  may  be  created  by  addin?  a  thin  layer  of  more  viscous  fluid  at  the 
inner  cylinder* 


Part  I*  Stability  equations  and  numerical  solution. 


He  use  cylindrical  coordinates  (r,6,z)  where  the  z~axis  is  the  axis  of  cylinders  of 


wl 


radii  R^  and  Rg.  He  consider  the  stability  of  a  circular  Couette  flow  of  two  fluids# 
lyin?  between  the  cylinders.  The  azimuthal  velocity  field  is  given  by  V^(r)  *  A^r  + 
Bj/r #  i  *  1f2,  where  i  *  1  refers  to  the  'inner*  fluid,  occupying  R^  <  r  <  D  and  i  -  2 
refers  to  the  'outer*  fluid#  occupying  D  <  r  <  R^*  The  unperturbed  interface  is  at 
radius  D.  The  angular  velocities  of  the  cylinder  are  the  viscosities  are  and 

the  densities  are  p4. 
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He  superi^ose  an  infinitesimal  disturbance 


(u(r)#v(r),w(r)ip<r))  exp(-i  ct  ♦  ies  ♦  in©) 


so  that  in  each  fluid#  the  tfavier-Stokes  equations  yield: 
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i(-  9  *  ^)v  ♦  2Au  -  -  i2E  +  v(l(rv'i'  -  (n2+1)2j  -  a2v  +  -2i|2)  (3) 

r  r 

vhtrt  v  -  u/P«  Incompressibility  yields 

(ru)'  ♦  inv  +  riaw  -  0  (4) 

The  int«r(«o«  Is  also  perturbed  and  tha  following  conditions  hold  at  r  -  D  (ssa  Nguyen  and 

Joseph#  1983).  [•]  danotas  tha  diffaranca  . 1  -  across  r  *  D. 

(1)  Continuity  of  velocity i  [u)  -  0,  [w]  •  0,  i(-  a  *  [vj  -  . 

D2 

(ii)  Continuity  of  itteinn 

iou  (Ml  +  [|»|^]  ■  0 

[p(f)’]  +  ^  tul  -  o 

D 

2 

1(-  a  ♦  2j)([pJ  -  2[uu'l)  +  utP'J  -  Tu(— — ♦  a2)  -  0 

D2 

2 

whars  T  is  tha  surface  tansion  and  P'  -  pV  •  Boundary  conditions  at  tha  solid  ares 

u  -  0#  v  -  0,  w  -  0  at  r  -  R-|  and  r2. 

We  now  dascriba  our  discretisation  schema.  Equations  (1)  and  (4)  are  used  to 
eliminate  w  and  p  from  aquations  (2)  and  (3).  A  Chebyschev  polynomial  expansion  (Orszag, 
1971)  is  used  for  u,v,w  and  p«  If  n+1  and  n  Chabyschav  polynomials  are  used  for  u  and  v 
respectively,  than  tha  total  number  of  unknowns  is  4n+2.  Equation  (2)  can  than  be 
truncated  after  tha  n-4th  degree  because  of  tha  presence  of  rVv  and  r3v* ' '  •  Equation 
(3)  should  be  truncated  after  tha  n-3rd  degree  because  of  the  presence  of  r3u'"  and 
r^v".  The  resulting  system  of  linear  equations  for  a  were  solved  with  an  DtSL  routine 
on  VAX/VMS  11/780  in  complex  double  precision*  The  computations  were  checked  against 
Table  2  in  Krueger,  Gross  and  DiPrina  (1966),  Hooper  and  Boyd's  asymptotics  for  large  n 
and  the  asymptotics  in  Part  II  for  large  a  . 

When  the  two  fluids  ere  identical,  the  presence  of  the  Interface  introduces  a 
neutrally  stable  eigenvalue  for  each  n  and  o  (called  the  'inter facial*  eigenvalue  by  Yih 
(1967))  in  addition  to  the  eigenvalues  for  one-fluid  flow  (called  'Taylor'  eigenvalues  in 
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Pert  III).  In  Part  III,  we  track  the  behavior  of  tha  intarfacial  and  Taylor  eigenvalues 


as  tha  viscosities,  densities,  aurfaca  tanalon  and  volume  ratio  are  changed. 


Part  II.  Asymptotic  analysis  for  abort  wave  diaturbancaa * 


Hooper  and  Boyd  (1983)  restrict  their  asymptotic  analysis  to  two  dimensions  with 


coordinates  (x,y)  where  x  is  the  direction  of  the  stream.  They  consider  disturbances  in 


armal  modes  proportional  to  exp  (iox)  for  large  a  or  short  waves ,  and  expand  the  stream 


function  and  the  interfacial  eigenvalue  in  powers  of  1/a  .  The  perturbation  problems 


which  arise  from  this  procedure  are  uniquely  solvable.  Since  our  stream  is  in  the 


asimuthal  direction,  we  must  replace  a  by  n/r.  The  results  of  Hooper  and  Boyd  apply 


when  n  is  large  and  centrifugal  effects  are  neglected.  Centrifugal  effects  considered 


here,  however,  play  the  same  role  as  gravity  in  their  analysis.  (We  note  that  a  factor 


a  should  multiply  the  gravity  term  below  (16)  in  their  paper).  For  large  n,  we  find 
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We  can  also  do  short-wave  asymptotics  for  disturbances  perpendicular  to  the  (r,0) 


plane.  We  consider  axisymmetric  (n  -  0)  disturbances  proportional  to  exp(iaz)  and 


introduce  the  following  dimensionless  variables:  R  ■  a(r-D),  T  -  2Bjt/D  ,  Z  -  bs/D. 

Assuming  now  that  the  disturbance  is  proportional  to  exp(ia1  CT  +  ia^Z)  where 
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this  transformation  ylalds  tha  following  six  conditions  it  R  *  0  as 
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Ths  aquations  of  motion  and  continuity  with  3/30  -  0  ara  now  expanded  about  R  ■  0 


for  larga  and,  as  in  tha  analysis  of  Booper  and  Boyd, 
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Wtor*  tha  seroth-ordar  velocity  satisfies  I>3u0  *  0  in  each  fluid.  Wa  find  that  c1  is 


determined  by  uq  where 
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of  tha  coafflciants  in  uq  can  ba  found  in  terms  of  tha  sixth  by  using  tha  intarfaca 


conditions  (1),  (2),  (4),  (5)  and  (6).  Condition  (3)  yialds  cQ  ■  0  and  an  aquation  for 
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Fro*  this  we  find  that 
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Part  III#  Numerical  roulti. 

Vi  compute  the  growth  and  decay  rates,  ia(o),  for  the  etabiXlty  of  Couette  flow  of 
two  f luida«  we  conaider  two  flow  reqiaes*  The  first ,  treated  under  (a)  and  (b)  below,  is 
flow  at  mil  Reynolds  nuabers.  Here,  if  either  fluid  filled  the  flow,  the  one-fluid  flow 
would  be  linearly  stable.  The  only  node  which  can  become  linearly  unstable  for  the  two- 
fluid  flow  is  the  interfacial  mode.  The  second,  treated  under  (c)  below,  is  flow  at 
hitter  Reynolds  nuabers  where,  if  the  outer  fluid  filled  the  flow,  the  one-fluid  flow 
would  be  at  a  critical  Taylor  nuaber  where  linear  stability  is  lost.  Here,  in  addition  to 
the  interfacial  eigenvalues,  the  eigenvalues  associated  with  the  one-fluid  flows  can 
become  unstable.  This  type  of  loss  of  stability  leads  to  bifurcation  and,  finally,  to  the 
tesselation  of  stable  (highly  viscous)  and  unstable  (Taylor  cells  in  the  low-viscosity 
liquid)  regions  observed  in  the  experiments  of  Joseph,  Nguyen  and  Beavers.  For  each  flow 
regime,  we  determine  which  arrangement  of  the  components  is  stable  and  the  volume  ratios 
of  the  stable  configurations. 

(a)  Stability  of  Couette  flow  of  two  fluids  for  low  Reynolds  nuabers  neglecting 
surface  tension  and  density  difference. 

We  compute  the  growth  rates  for  the  following  range  of  variables i 
0,  -  1,  Q2-  °,  M2-  a-  0.2  to  6,  R1  -  1,  R2  -  2,  p,  «  p2  -  **  the  Reynolds  nuaber 

Re  -  V(R1)  (R2  -  R1)/v1  ranges  from  0.5  to  5,  a  ranges  from  0.01  to  50  and  n  from  0  to 
50.  Under  these  conditions,  we  find  that  the  configuration  with  a  sufficiently  thin  layer 
of  the  less  viscous  fluid,  situated  next  to  either  cylinder,  is  stable. 

The  response  to  long  waves  (small  a  and  low  n)  is  as  follows.  The  axlsyametric 
mode  becomes  insignificant  as  o  ♦  0  since  in  that  limit,  there  is  no  disturbance.  For 


a  <  0.1  and  small  Reynolds  numbers#  the  growth  rate  im<o)  is  proportional  to  a  Re  when 
n  *  0  and  to  Re  when  n  #  0  .  The  growth  rates  in  this  asymptotic  range  are  shown  in 
figures  1  and  2.  In  figure  1,  the  leas  viscous  fluid  is  situated  next  to  the  inner 
cylinder  and  hence  the  modes  displayed  are  stable  (im  (o)  <  0)  if  the  interface  is  close 
enough  to  Rj  *  1.  The  situation  is  reversed  in  figure  2*  In  both  figures#  modes  10  and 
20  show  the  short-wave  asymptotic  behavior  in  which  the  stable  range  of  interface 
positions#  as  well  as  the  maximum  growth  rates#  diminish  with  n  (or  a  )•  Disturbances  of 
the  stable  configurations  with  the  less  viscous  fluid  inside  ha  juo larger  decay  rates 
than  that  with  less  viscous  fluid  outside. 

Trends  similar  to  those  exhibited  in  figures  1  and  2  are  *  >ited  in  figures  3  and  4 
for  a  -  1.0.  The  stable  range  of  Interface  radii  is  slightly,  ,»ot  greatly, 
reduced.  Figure  4  clearly  shows  that#  for  m  -  2#  the  dependence  of  im(o)  on  n  at  modes 
20  and  40  scales  with  1/n2  over  most  of  the  Interface  positions.  In  both  figures#  the 
relative  errors  of  the  asymptotic  values  at  D  -  1.5  fall  from  about  50%  at  mode  9  to  8%  at 


Figures  5  and  6  give  growth  rates  for  a  ■  10.  For  fixed  small  values  of  n  and  large 


a#  we  find  that 


13  (1-a)2  2.,  D2 

c.  ~ - m  (1 - =■)  and 

0  <H*)Z  Rj 

C1  1  -4 

im(o)  -  (— j)  +  0(a  *)  . 

D  W1  a 


In  figures  5  and  6#  modes  n  ■  1  to  3  lie  in  between  modes  0  and  4  and  the  growth 
rates  of  all  the  modes  between  0  and  4  are  numerically  close.  These  figures  display  some 
qualitative  features  of  large  a  asymptotics#  but  a  *  10  is  not  high  enough  to  be  in 


the  short-wave  asys$totic  range  for  m  *  0.4  or  2.  in  addition,  the  lar  the  viscosity 

difference  m#  the  lower  is  the  value  of  o  at  which  this  asymptotic  range  is  attained. 
For  example#  at  m  -  6#  the  relative  errors  range  from  30%  at  a  ■  10#  16%  at  a  *  20 
and  8%  at  a  *  40#  whereas  at  m  -  2  (figure  6),  these  errors  are  doubled  and  at  m  *  0.4 
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(figure  5),  they  are  quadrupled.  As  noted  by  Hooper  and  Boyd,  the  short-wave  asymptotics 
breaks  down  when  the  interface  is  too  close  to  R-j  or  R2. 


From  our  computations,  we  conclude  that  the  largest  growth  and  decay  rates  arise  in 
the  order-one  range  of  a  for  medium  n  and  for  small  m.  For  example,  mode  9  at  o  -  1, 
m  -  0.2,  attains  triple  the  growth  rate  attained  at  m  -  0.4  (figure  3)  and,  in  turn,  that 
mode  at  m  *  0.4  attains  a  larger  magnitude  than  at  m  ■  2  (figure  4).  We  may  also  conclude 
from  a  coaparison  of  decay  rates  that  the  stable  flows  with  thin  fluid  inside  are  *more’ 
stable  than  those  with  thin  fluid  outside* 

(b)  Stability  of  Couette  flow  of  two  fluids  for  low  Reynolds  numbers.  The  influence 
of  surface  tension  and  density  differences. 

Surface  tension  stabilises  short  wave  Interfacial  disturbances  and  destabilises 
longer  waves.  Centrifugal  forces,  in  the  absence  of  surface  tension,  will  produce 
stability  if  the  more  dense  fluid  is  outside.  However,  with  surface  tension,  it  is 
possible  to  achieve  stability  when  the  denser  fluid  is  inside.  This  can,  of  course,  only 
happen  if  the  centrifugal  force  is  not  too  large  and  gravity  is  neglected.  Under  these 
conditions,  if  surface  tension  is  large  enough  to  stabilise  the  short  waves  but  not  so 
large  that  the  long  waves  are  unstable,  then  stability  is  possible  at  all  a  and  n  with 

the  denser  fluid  inside.  One  exa^>le  is  T  *  1,  r  *  2,  m  -  2,  D  ■  1.9,  u2  “  1*  P2  *  1  • 

Figure  7  shows  a  graph  of  -im(a)  versus  a,  showing  stability.  Figure  8  shows  a  graph 
of  im(0)  versus  a  at  sero  surface  tension,  showing  that  modes  become  unstable  for  large 
a  . 

( c )  Stability  of  Couette  flow  of  two  fluids  near  a  critical  Taylor  number.  Zero 
surface  tension  and  density  difference* 

We  study  the  onset  of  Taylor  insteVMity  associated  with  modes  n  *  0  and  1.  We  also 
study  the  higher  modes  (n  >  1)  for  which  the  Taylor  modes  are  stable.  In  the  classlclal 
Taylor  problem  for  one  fluid,  the  most  unstable  mode  is  the  axisymmetric  (n  ■  0)  one. 

4 

The  Taylor  number  is  defined  as  T^  •  dAft^R^  -  )  /v  where 

A  -  (fl.-R  2  "  0  R  2)/(r  2  -  R  2).  When  fl2  -  0,  fl.  -  1,  R1  -  1,  R,  -  2,  p.  -  p3  -  1,  the 


1549.59  and 


2 

a  xi  symmetric  Mode  becoaei  unstable  for  Ta  close  to  Tcr  where  Tcr(R1/R2>  ■ 

e 

a  ■  3.16  (Krueger,  Gross  and  DiPrima,  1966).  We  fix  v2  *  0.0146  and  vary  so  that  if 
the  outer  fluid  occupied  the  whole  annulus,  (l.e.  D  ■  R^ )  the  growth  rate  for  the 
axisyssMtric  Taylor  Mode  is  near  criticality.  We  note  that  this  choice  of  implies 

that  the  i*(o)  for  the  axisyametric  Mode  passes  through  zero  when  D  ■  as  in  figures 
9-12.  With  this  choice  of  paraneters,  we  study  how  the  growth  rates  vary  with  D. 

We  consider  two  situations t  the  less  viscous  fluid  is  at  the  inner  (m  <  1,  figures  9 
and  10)  or  at  the  outer  (m  >  1,  figures  11  and  12)  cylinder.  Intuition  would  suggest  that 
when  m  >  1,  the  flow  will  becoae  More  and  More  stable  as  D  increases  because  of  the 
presence  of  increasingly  larger  amounts  of  More  viscous  (stable)  fluid.  This  expectation 
is  not  realized.  Figures  11  and  12  show  that  various  Modes  become  unstable  as  thick  fluid 
is  added.  Similarly,  intuition  would  suggest  that  when  m  <  1,  we  should  have  Instability 
for  increasing  D  because  more  and  More  thin  fluid  replaces  thick  fluid.  Figures  9  and  10 
show  that  we  actually  stabilize  the  flow  by  adding  less  viscous  fluid  near  the  inner 
wall.  This  stabilization  near  D  -  1  is  associated  with  the  stability  of  narrow  layers 
near  the  inner  cylinder  and  could  be  called  'lubrication*  stabilization  associated  with 
the  layer  of  thin  fluid  on  the  inner  cylinder. 

A  new  feature  close  to  or  above  a  critical  Taylor  number  is  that  the  im(o)  for  the 
interfacial  eigenvalue  need  not  be  single-valued.  That  is,  the  graph  of  im(o )  versus  D 
for  an  interfacial  eigenvalue  which  begins  at  D  ■  R^  with  im(o)  *  0  can  proceed  to  Match 
to  a  Taylor  eigenvalue  at  D  -  R2  and  a  second  branch  satisfying  im(o)  ■  0  at  D  *  Rj  can 
match  to  a  Taylor  eigenvalue  at  O  *  R^.  Figures  10  and  11  show  mode  1  to  have  such 
branches.  In  figure  10,  the  Taylor  eigenvalue  for  mode  1  is  unstable  at  D  *  Rj  and  in 
figure  11,  it  is  stable  at  D  *  R2*  The  behavior  of  the  higher  (n  >  1)  Modes,  for  which 
the  Taylor  modes  are  very  stable,  is  as  described  in  part  (a). 

For  n  ■  0,  the  equations  yield  a  real-valued  problem  for  io,  ip,  iu,  iv  and  iw. 

Hence,  io  is  either  a  real  number  or  appears  in  complex  conjugate  pairs.  In  the  latter 
case,  the  two  eigenvalues  have  equal  imaginary  parts.  This  behavior  is  shown  in  figures  9 
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to  12*  For  axaapla,  in  figura  9f  for  Rj  <  D  <  1.3,  tha  n  »  0  aiganvaluas  ara  In 
con jugata  pairs.  Naar  D  *  1.3,  tha  in(o)  splits.  Ona  branch  bacoaas  incraasingly 
unstabla  and  at  D  ■  Rj  is  tha  unstabla  Taylor  aiganvalua  for  ona-fluld  flow  with 
v  -  0.0132.  Tha  sacond  branch  bacomas  stabla  for  1.3  <  D  <  ^  and  *inca 
i*(o)  ■  0  at  0  ■  Rj  this  branch  is  tha  intarfadal  aiganvalua. 
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ALPHA=0.1 


Figure  1 


Growth  rate  versus  Interface  position  with  azimuthal  wave  number  (n)  as  a  parameter* 

T  •  0,  r  If  ■  ■  0.4;  The  less  viscous  fluid  is  on  the  inner  cylinder*  Negative 
im  (o)  corresponds  to  stability*  When  surface  tension  is  absent/  the  flow  is  unstable 
at  any  D  (  ^  or  Rj  )  if  n  is  large  enough*  Mode  0  is  insignificant  under  graph 

scales* 
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ALPHA-0.1 


INTERFACE  RADIUS  (D) 


Figure  2 


Growth  rat*  versus  interface  position  whan  tho  loss  viscous  fluid  is  on  ths  outsr 
cylinder.  T-0,  r  ■  1,  a  ■  2.  Ths  stsbls  sodas  nasr  ths  outsr  cylinder  have  lass 
stability  than  ths  stsbls  nodes  near  ths  inner  cylinder  (ef  figure  1)  because  ths  decay 
rates  of  stsbls  disturbances  are  such  snaller. 
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Figure  3 


Growth  rate  curves  when  the  less  viscous  fluid  is  inside*  T  *  0,  r  ■  1#  ■  -  0,4* 


mi 


WV. 


Figure  4 

Growth  rate  curves  when  the  less  viscous  fluid  is  outside*  T-0,  r»1,  m-2*  The 
decay  rates  of  stable  disturbances  are  an  order  of  Magnitude  smaller  than  in  figure  3. 
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Figure  9 


Qrowth  r**M  lot  various  aodee  when  the  parameters  ere  close  to  critical  for  Taylor 
Instability.  T  »  0,  r  »  1,  s  ■  0.9.  Thick  fluid  lies  next  to  the  outer  cylinder.  The 
addition  of  thin  fluid  at  the  inner  cylinder  surprisingly  stabilises  the  flow. 


0.04 


ALPHA=3.16 


Figure  11 


Growth  rates  when  the  less  viscous  fluid  lies  next  to  the  outer  cylinder.  T  -  0,  r  ■  1 , 

■  ■  1.08.  The  amount  of  less  viscous  fluid  decreases  as  D  increases  but  various  modes  are 
unstable  except  when  most  of  the  gap  is  occupied  by  the  more  viscous  fluid.  At  D  -  Rj, 
mode  0  is  at  a  critical  Taylor  number  and  mode  1  is  slightly  below. 
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